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Abstract 

Taylor  and  Yates  have  measured  atmospheric  absorption  in  the  infrared 
over  long  horizontal  paths.  They  published  thirteen  spectra  which  covered  four 
path  lengths  from  0. 3  to  27. 7  kilometers  and  seven  water  vapor  concentrations 
from  0. 10  to  2. 3  pr  cm/km.  This  report  describes  a  method  to  fit  tiveir  spectra 
to  a  model  in  the  3-  to  5-micrometer  interval.  With  this  method,  particulate 
scattering  coefficients,  band  parameters  of  the  Goody  model  for  H20  absorption 
and  band  parameters  for  the  combined  absorption  of  all  other  molecular  apecies 
were  extracted. 
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I.  Introduction 

Considerable  effort  has  beer,  expended  over  the  past  several  years  to 
provide  accurate  atmospheric  transmission  models  which  are  analytically  trac¬ 
table  and,  thereby,  readily  adaptable  to  digital  computers.  Anding  [1]  has  pro¬ 
vided  a  discussion  of  present  computational  models.  The  well-lnown  experi¬ 
mental  data  of  Taylor  and  \ates  [2,3]  are  widely  used  to  compare  various 
models.  Since  the  models  often  disagree  considerably  in  particular  regions  of 
the  spectrum,  it  was  felt  uesirable  to  fit  the  data  of  Taylor  and  Yates  to  a  band 
model  n  a  spectral  interval  of  immediate  interest  to  several  work,  '  groups. 

This  report  discusses  a  method  to  derive  a  set  of  band  parameters  from 
the  measured  transmission  data  of  Taylor  and  Yates  in  the  3-  to  5-micrometer 
spectral  interval.  Taylor  and  Yates  measured  the  extinction  of  light  which  had 
traversed  various  distances  from  a  carbon  arc  t>earchlamp  to  a  monochromatic 
receiver  over  land,  across  Chesapeake  Bay,  over  water  parallel  to  the  bay 
shoreline,  and  between  two  Hawaiian  mountain  tops  at  10,  000  feet. 

Two  effects  cause  attenuation  of  the  light  —  scattering  by  small  water 
droplets  and  other  particles,  and  absorption  by  certain  molecular  species,  such 
as  H20,  C02,  NzO,  NH4,  and  CO.  The  overall  transmission  for  a  given  wave¬ 
length  is 


T<\)  =  T  (\)T  (X)  , 

a  “ 

where  a  and  s  refer  to  absorption  and  scattering,  respectively. 

Since  the  amount  of  aerosol  was  not  recorded  and  the  aerosol  cross 
sections  vary  with  meteorological  conditions,  to  separate  T  and  T  would  be 

E  S 

impossible  were  it  not  for  spectral  differences  between  the  two.  Lines  accent 
the  T  spectrum,  but  T  is  only  weakly  wavelength  dependent.  Sharp  variations 

3  S 

j  ■  the  spectrum  can  be  attributed  to  absorption,  and  overall  attenuation  to 
scattering.  Because  the  absorption  by  many  overlapping  spectral  lines  can  be 
mistaken  for  scattering,  the  separation  of  T  and  T  will  be  slightly  ambiguous. 

Since  two  parameters  are  varied  in  the  Taylor  and  Yates  data  —  the 
amount  of  water  vapor  in  the  path  and  the  length  of  the  path  — .  one  can  expect 
to  separate  the  molecular  absorption  into  that  of  water  vapor  and  that  of  all 
constituents  of  constant  concentration.  Unavoidably,  the  absorption  by  C02, 
N20,  CH4,  and  CO  is  combined,  the  sum  of  which  will  be  attributed  to  a  single 
species,  X02.  The  amount  of  X02  was  chosen  to  be  32  atm  cm/km,  which  is 
the  average  concentration  of  C02. 
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The  Goody  model  [4]  was  selected  for  its  computational  simplicity.  It 
is  described  by  the  expression 

1/ 

T  =  exp{-(S/i)i  (1  +  2(S/27ra)Pi  )/2}  . 

cl 


In  this  model  the  spectral  lines  are  assumed  to  be  located  randomly  in  frequency 

with  a  mean  spacing  d,  half  widths  a,  and  strengths  exponentially  distributed 

with  a  mean  intensity  S.  The  amount  of  absorbing  material  is  t ,  and  P  is  the 

total  pressure  in  atmospheres.  Temperature  effects,  which  would  introduce 

extra  parameters  into  the  model,  were  not  considered.  This  model  is  physically 

applicable  for  molecules  with  low  symmetry.  A  modei  with  a  uniform  line 

spacing  describes  linear  molecules  such  as  C02  and  N20  more  accurately.  All 

models  predict  T  cc  i  for  small  *  and  T  a  #  for  large  i  .  If  the  data  are 
R  R 

limited  to  one  of  these  asymptotic  regions,  the  curves  of  growth  can  be  described 
by  the  Goody  model  as  well  as  any  physically  accurate  model.  The  band  param¬ 
eters  are  not  uniformly  accurate  over  the  entire  band.  In  regions  of  weak  H20 
absorption,  it  is  possible  to  determine  the  X02  parameters  precisely,  but  when 
H20  absorption  masks  X02  absorption,  to  deduce  X02  parameters,  albeit  roughly, 
is  difficult  even  though  its  absorption  may  be  as  strong  as  in  the  former  case. 


The  goal  of  this  empirical  approach  is  to  separate  the  transmission  into 

three  components,  T  -  T  T„  _T„  .  The  method  to  extract  T  is  discussed  in 

s  H20  X02  a 

section  2. 


The  algorithm  to  separate  T  from  T  is  described  in  section  3. 

The  least  information  that  this  empirical  approach  yields  is  a  set  of  parameters 

that  reproduce  the  Taylor  and  Yates  spectra.  Reproducing  the  fitted  data  is  not 

difficult  since  ao  error  in  T  can  be  corrected  accidentally  by  an  opposite 

H2U 


error  in  T, 


At  best  the  method  can  determine  T 


th2o’  and  txo2 


X02  s' 

separately.  If  known  N20  absorption  appears  as  an  increase  in  the  X02  param¬ 
eters  and  not  in  the  H20  parameters,  the  method  is  successful. 


2.  Particulate  Scattering 

A  scattering  coefficient  k,  defined  by  T  =  exp  (-k£)  where  I  is  the 

s 

path  length,  was  found  by  the  procedure  described  below.  If  a  wavelength  \  for 

which  there  was  no  molecular  absorption  existed,  then  the  scattering  could  be 
found  by  T  =  t/T  =  T(\0)/l=  T(A0) .  For  the  short  path  lengths  this  proce- 

dure  is  possible.  However,  T  *  1  for  all  wavelengths  of  the  Taylor  and  Yates 

spectra  over  the  longer  path  lengths,  and  a  different  procedure  must  be  used. 
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ft 


From  the  fine  resolution  spectra  of  Howard  [5] ,  one  finds  only  extremely  weak 
lines  in  the  interval  3.96  to  4.01  micrometers.  In  this  interval,  the  transmis¬ 


sion,  which  is  attenuated  by  the  wings  of  C02,  has  the  form  T  =  exp^-aPf 


xo2. 


a  V 

The  coefficient  a  was  derived  from  spectra  of  short  path  lengths  for  which 

T  =  T^).  The  unknown  scattering  was  found  by  T  =  t/t  =  Texp/aPiv 
e  say  aU 

This  method  relies  on  a  single  spectral  point  tn  derive  the  scattering  for  the 

entire  spectrum;  a  global  procedure  is  more  reliable.  A  band  absorption  model 


>> 


was  calculated  with  the  set 


T  ' 
s 


fT 

I  •/ 


A  new  set 


was  calculated  by 


where  i  indexes  the  spectral  points  and  T 


■CO 


was 


calculated  by  the  model.  The  final  band  absorption  model  was  calculated  with 


Table  I  shows 


|Tg|  and  Since  the  absorption  is  temperature 


dependent,  neglecting  temperature  effects  in  the  model  would  introduce  some 

scatter  in  the  values  of  T  .  An  attempt  to  improve  the  sums  of  errors  of  the 

s 

fit  by  allowing  a  wavelength  dependent  k  yielded  slopes  in  the  interval 
±0. 02  k/ micrometer,  but  failed  to  improve  the  error  enough  to  indicate  the 
slopes  were  real. 


•TABLE  I.  SCATTERING  COEFFICIENTS 


Path  Length  (km) 

H20  (pr  cm) 

T 

8 

T  » 
s 

k  (km~‘) 

0.3 

0.1 

1.00* 

1.00* 

0.3 

0.6 

1.00* 

1.00* 

5.5 

1.4 

0.79 

0.79 

WmmBm 

5.5 

2.2 

0,75 

0.  79 

5.5 

4.2 

1.0C 

0.97 

0.005 

5.5 

9.4 

0.97 

0.93 

0.  014 

16 

5.2 

0.58 

0.63 

0.028 

16 

6.7 

0.  67 

0.67 

0  325 

16 

15 

0.85 

0.  84 

0.011 

16 

28 

0.63 

0.  66 

0.  025 

16 

37 

0.95 

0.91 

0.  006 

27 

10 

0.60 

0.55 

0.  022 

27 

20 

0.69 

0.70 

0.013 

*By  definition. 


3,  Molecular  Absorption 


The  Goody  model  can  be  rewritten 


In  regions  of  no  X02  absorption,  c  and  d  should  be  zero.  *  However,  a  fit  would 
likely  give  spurious  non-zero  c  and  d  coefficients.  Contractions  of  the  model 
must  be  considered.  Let  (a,  b,  c,  0)  denote  the  model 


The  notation  (a,  b,  0,  0) ,  etc. ,  is  defined  similarly. 


An  algorithm  was  devised  to  make  a  fit  with  the  optimal  number  of 
coefficients.  Hopefully,  the  algorithm  could  decide  to  use  (a,  b,  0,  0)  when 
X02  absorption  was  negligible.  First,  the  beet  one-coefficient  model  was 
chosen  from  the  models  (a,  0,  0,  0),  (0,  b,  0,  0),  (0,  0,  c,  0),  and  (0,  0, 

0,  d).  Suppose  (0,  b,  0,  0)  had  the  least  sum  of  squares  of  errors.  The  best 
two-coefficient  model  was  chosen  from  all  possible  combinations  of  two 
coefficients.  Assume  it  was  (a,  0,  c,  0).  From  (a,  b,  c,  0)  and  (a,  0,  c,  d) , 
the  best  three-coefficient  model  was  chosen,  say  (a,  0,  c,  d) .  Intuitive  criteria 
picked  the  optimal  model  from  the  three  best  models;  in  this  example.  (0,  b, 

0;  <t),  (a,  0,  c,  0),  and  (a,  0,  c,  d) .  Since  the  three-coefficient  model  would 
mirror  the  data  most  precisely  —  fortuitously  in  most  cases  —  the  least  errors 
condition  is  insufficient  to  reject  accidental  fits.  The  Goody  model  precludes 


1/ 

negative  coefficients  because  -log  T  =  (  Pi  /a)  2  when  i  is  large  and 
1/  a  '  ' 

-log  T  =  (i2/b)  2  wheni  is  small.  The  first  test  utilizes  this  fact.  Criterion  1 


is  to  rejc>.'  models  with  negative  coefficients.  However,  spurious  models  with 
positive  coefficients  would  pass  criterion  1.  In  most  oases  the  fortuitous  mode! 
improves  the  error  only  slightly  over  the  true  model  of  fewer  coefficients. 
Criterion  2,  to  minimize  Z(errors)2/ (No.  of  data  points  —  No.  of  coefficients) , 
would  veed  out  the  unwanted  models.  Perhaps  statistical  theory  requires  some 
other  factor,  but  this  question  was  not  explored. 


*The  number  l/ 0  is  interoreteri  to  be  n  so  that,  o  =  d  =  0  instead  of 
c  =  d  =  oo  indicates  absence  of  X02. 
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In  the  C02  wings  the  model  was  replaced  by 

-“*Ta=  P<X02/^  '  (2) 

The  pressure  factor  is  P  in  the  Lorentz  theory  of  line  profile,  and  Winters, 
Silverman,  and  Benedict  [6]  have  confirmed  it  for  the  far  wings  of  the 
4. 3-micrometer  C02  band. 

Since  the  distance  between  the  source  and  detector  was  measured,  the 
perfect  gas  correction  for  temperature  and  pressure  was  made  to  reduce  i 

XU  2 

amounts  to  a  standard  sea  level  reference,  by  using 

iX02(T,P/iX02(T°’Po)  =  (p/po)(T</T)  . 

t 

All  fits  were  done  by  the  least  squares  method,  to  minimize  errors  in 
-log  T  .  Datum  was  not  used  if  T  s  o.  98  or  T  ^  0.  05.  The  average  error 

S  2L  <1 

was  <  f  AT  /T  \  2  >  =  0.  02. 

\  i  v 

4-.  Calculated  Band  Parameters 

The  band  parameters,  which  were  calculated  with  the  scattering 
transmission  set  /t  '( ,  are  listed  in  Table  II.  The  wavelengths,  which  are 

l  8  ) 

meant  to  identify  spectral  features  such  as  maxima  and  minima,  have  errors 
no  larger  than  0.  005  micron.  The  model  is  capable  of  reproducing  the  Taylor 
and  Yates  spectra  sufficiently.  Figure  1  compares  the  fit  the  model  generated 
and  one  Taylor  and  Yates  spectrum.  These  coefficients  should  not  be  extrap¬ 
olated  beyond  the  limits  (0.  3  km,  17  km)  ac  sea  level  and  (0. 1  pr  cm  H20, 

37  pr  cm  H20) . 

To  decide  if  the  algorithm  separated  X02  from  H20,  one  transmission 

curve  was  compared  to  the  laboratory  data  of  Howard,  Burch,  and 

Williams  [7,8]  for  individual  species.  Since  the  conditions  of  their  artificial 

atmosphere  did  not  match  those  of  the  Taylor  and  Yates  data,  comparisons 

should  decide  qualitatively  how  well  X02  was  separated  from  H20,  and  precise 

agreement  is  not  expected.  In  Figure  2,  -log  T  was  plotted  for 

H2vj 

_  =  9.4  pr  cm,  and  -logT  for  1  =  176  atm  cm.  The  latter,  the 

112^  AU2  XU2 

lower  curve,  was  shifted  down  by  a  factor  of  10  to  offset  the  two  curves.  Super¬ 
imposed  on  this  graph  are  several  transmission  curves  of  Howard,  Ejrch,  and 
Williams,  i  rom  3.0  to  3.5  micrometers.  H-O  is  the  H.  absorber.  In 
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the  2^2  band  at  3.  2  micrometers,  the  CH4  band  center  barely  appears  against 
strong  H20  absorption,  but  the  baud  edge  is  seen  clearly  where  H20  absorption 
is  weaker.  From  3.  5  to  3.  8  micrometers  the  HDO  band  is  a  prominent  feature 
in  the  H20  curve.  Most  HDO  peaks  were  matched  with  HDO  line  clusters  of 
Howard  [5].  The  X02  absorption,  in  this  interval,  may  be  attributed  to  CH4 
lines.  From  3. 8  to  4. 2  micrometers  there  are  two  N20  bands  and  the  feature¬ 
less  C02  wings.  The  N20  and  C02  bands  dominate  from  4. 36* to  4.6  microns. 
From  4.  6  to  5.  0  microns  most  of  the  H20  peaks  were  identified  with  clusters  of 
H20  lines  in  the  solar  spectrum  of  Shaw,  Chapman,  Howard,  and  Oxholm  [9], 
Half  of  the  4. 8-micrometer  C02  band  was  masked  by  stronger  H20  absorption. 
The  CO  band  was  not  resolved  clearly. 

A  computer  program  that  computes  transmission  spectra  with  the  band 
model  of  Table  n  is  listed  in  the  appendix. 


TABLE  II.  BAND  PARAMETERS 


y 

d'2  (atm  cm) 

*/«  '/  v  y 

"(pr  cm-atm)  2  *g  2  (atm  cm-atm)  c^J 


Wavelength 

/ 

(|im) 

b/J  (pr  cm) 

3.000 

3.010 

7.415F-01 

3.020 

3.030 

3.040 

3.050 

2.624E-01 

3.060  * 

3.741 E-C 1 

3.070 

8.916F-01 

3.080 

3.090 

7 . 058F-0 1 

3.100 

7.485E-01 

3.110 

..2.451E-01 

3.  120 

6.864E-01 

3.130 

8.164Z-01 

3.140 

3.  160 

3.497E-00 

3.170 

2 • 232E-00 

3.200 

1. 1C6E-01 

3. 216  ‘ 

2 . 199E-01 

3.223 

1 .587F-01 

3.230 

4.634E-01 

3.240 

3.250 

fH 

o 

1 

LU 

o 

r-f 

• 

3.260 

3.270 . 

7 • 41 3E-00 

3.290 

3.300 

5.687E-01 

3.310 

3.327 

1.140E-00 

3.333 

'  3.340“ 

9. 123E-01 

3.350 

8.507E-01 

3.370 

3.  380 

3.400 

2.281F+01 

3.402 

3^420 

3.427 

3.433 

3.440 

3.  177E  +  01 

3.445 

3.450 

3 . 097E+0 1 

3.460 . 

4.831E+01 

3.470 

4.243E+01 

3.490 

4. 250EJ-01 

3.500 

3.867E+01 

3.514 

6.068E+01 

6.936E-01 
7 . 694£~oi 
8.273E-01 
9.229F-01 
7.297E-01 
5. 133  E-01 
6.699E-01 

8.591E-01 


4.975E-01 
8.095E-C1 
9.472E-01 
8.712E-01 
1  .464E-00 
1  .412E-00 
3.825E-01 
4.614F-01 
4.460E-01 
8 . 925E-01 
7.676F-01 
1 •  105E-0G 

9.715E-01  3.4C6E+02 

7 . 172E-G1 
7.62IF-01 
8.914E-01 
1.519E-00 
1 . 166E-00 
1.162^-00 
1 •022E-00 
i .818E-00 
1 • 866E-00 

2 . 526E-QO 
5  »  72 7E-00 
5.092E-00 
5.083E-00 

5 .906F-00 


y. 

(atm  cm-atm)  ‘ 


4.495F+01 


1.602E+01 


6.063E+01 
.  • 1 87E+01 


8.221E+01 

9.850E+01 
1 .362E+02 
1.311E+02 
2.695E+02 
2 • 199E+02 
2 • 344E+02 
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TABLE  II.  BAND  PARAMETERS  (Continued) 


Wavelength 

(fxm) 

3.525 

3.535 

_3_«  54Q _ 

3.550 
.555 
.  560 
.570 
.580 

*  585 _ 

.600 
.515 
.625 
3.635 
3.645 

3.650 _ 

3.655 

3.670 

3.700 

3.710 

3.720 

3.730 

3.745 

3.755 

3,770 

3.790 _ 

3.800 

3.810 

3.820 

3.830 

3.840 


y 

cr2  (atm  cm) 

1/  1/  y  y 

b/2  (pr  cm)  a/2(pr  cm-atm) ' 5  (atm  cm-i 


y,  1/ 

atm)  c  z  (atm  cm-atm)  2 


3.870 

1. 

3.880 

8. 

3.910 

3.930 

6. 

3.940 

1. 

3.950 

1. 

3.960* 

4.000* 

4C922E+01 
4. 382E+01 
1 ,5266+01 
1.737E+01 
1. 2036  +  0 1 
2.257E+01 
1 . 235E+0 1 
2.420E+01 
2.557E+01 
5.271E+01 
3.396E+01 
7 . 459E+0 1 
_4_.038£+01 
9.098E+0 1 
_4.293E+0L 
4.945E+01 
3.250E-00 
9<»  463E+01 
_3.471E+01 
7. 476E+01 
_3. 1116+01 
6.232E+01 
2.450E+01 
4.970E+01 

77233E+02 

1.535E+01 
8,  637E+01 
3.307E+01 
1.051E+02 
4.636E+01 


2 . 529E+03 
3 . 068E+03 


3.526E-00 

3.391E-00 

3.470E-00 

6.435E-00 

4.347E-00 

6.051E-00 


1.768E+02 

2.777E+02 

4.338E+02 
6.263E+02 
1 . 832E+02 
1  •  533E+02 
2 . 563E+02 

w  •  J  O  i  ' 

3. 198E+02 


2.409E-00 _ 


5. 850E-00 
3._932E-Q0 


3.727E+01 


4.050# 

4.060* 

4. 100* 

4. 140* 

4.150* 

4,160* 

4.180 

4.190 

4.350 


7.1 10E+03 


3.424E+03 
3.717E+03 
2.11 2E+03 
1.681E+03 
7.358E+02 
2.646E+03 
1 . 152E+03 
1.397E+03 
_ 1.092E+03 
2 • 069E+03 
.  1.413E+03 
6.200E+02 
-±m&93E.±QZ. 
3 • 562E+02 
2.139E+02 
K866E+02 
1.798E+02 
9.545E+01 


3.79.9E+02 


4.750E+02 


4.816E+02 


2.965E+01 


7.772E+01 

8.354E+01 


♦Model  of  Equation  (2) ,  otherwise  model  of  Equation  (1) . 
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TABLE  II.  BAND  PARAMETERS  (Concluded) 


Wavelength 

(/on) 

4.377 

4.380 

4.450 


U 

(  i  d/J  (atm  cm) 

b/j  (pr  cm)  a'*(pr  cm-atm)  ^  *g'/*  (atm  cm-atm)  (atm  cm-atm)  ^ 


6.101E  00 
5.356E  00 
1.499E  00 


4.490 

2 . 567E+0 1 

1 . 173E  +  02 

4.500 

1 . 285E+0 1 

7.415E+01 

■“47520 

1.060E+01 

8. 1 15E+01 

4,530 

6.491E-00 

9.816E+01 

4.  554' 

3.266E-C01 

9.1Q1E+01 

8.626E-00 

4.  560 

2.577E+01 

7.249Ef01 

9.066E-00 

4.587 

2.544E+01 

3.091E+01 

4.595 

2 • 266E+0 ] 

2.676E+01 

4.605 

2 . 397E+0 1 

5  •  136E+01 

4.615 

.  2.113F-00 

4.630 

1  *  496E  +  0 1 

1.173E+02 

4.637 

1.378E+01 

7.234E+01 

4.643 

1  • 355E+0 1 

1.490E+02 

4.650 

2.121E-00 

2.041E-00 

•  4.660 

i.359E+01 

4.670 

1.406E-00 

4.660 

1 • 684E+0 1 

4.  690 

1.^77E+01 

2  «  027E+02 

4. '7  00 

1.343E-00 

1.834E-00 

4.715 

_ 1 «  348E  +  0 1  ... 

1.672E+02 

4,720 

1 .633E-00 

4.730 

1  • 363E+01 

4,  740 

1.537E-00 

2.094E-00 

4.745 

1 .  108E+01 

1 . 625E-00 

1 • 984E+02 

' 4.755 

3. 596E-00 

1 .630E-00 

4.784E+02 

4.760 

_ 3.270E-00 

1.7UE-00 

4.780 

1 -032E-00 

4.790 

9.529E-00 

1 . 447E+02 

4.800 

6, 365E-00 

1 • 385E+01 

4.81G 

4.715E-00 

3.338E+01 

'4.820 

8.377E-00 

4. 189E+01 

4.830 

6.575E-01 

5TB  415 

- 7;~575E~01'  " 

9.834E-01 

1  *  022E+02 

4.855 

2.674E-00 

1.470E-00 

4.860 

2  v  78  IE-00 

1.306E-00 

4.870 

2.909E--00 

1.319E-00 

4.890 

6.757E-01 

4.  Q00 

4, 392E-00 

4.920 

7. 77 7 K -00 

4.930 

1.116E-00 

4.940 

7.341E-01 

7.32iE-01 

4.955 

1 .871E-01 

4.965 

2.903E-00 

4.975 

7.618E-01 

4.990 

’  Zi  .  PAftF-Orv 

5.000 

3.638E-01 

7.453E-01 
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Appendix 

Fortran  IV  Program  to  Calculate  Tran*rni»sion  Spectra 


The  Fortran  IV  program  listed  here  computes  transmission  spectra 

with  the  band  parameters  :>f  Treble  II.  The  necessary  input  data  are:  a)  the 

amount  of  water  in  the  path  in  precipitable  centimeters,  b)  the  sea  level  path 

length  in  kilometers,  c)  the  pressure  in  atmospheres,  and  d)  the  amount  of 

scattering.  A  sample  set  of  input  and  output  is  shown  for  i  =  34.5  pr  cm, 

H2O 

path  length  15  km,  1  atm  pressure  and  Tg  =  1.  To  calculate  several  spectra 

at  a  time,  additional  input  cards  are  inserted,  one  for  each  spectrum.  This 
program  was  written  for  the  IBM  7094  computer. 


SIBFTC  HCLEB 

_ DIMENSION  W( 144) .C0EF(4,144) ,T(  144) ,L(4) tA(4I . FLAG M 4& 1 »F 11 44  1 _ 

DIMENSION  WAVE(48,3),TAU(48,3) 

_ LNT.E.<ld.fLf _ 

REAL  L 

_ LOGICAL  FLAG _ 

EQUIVALENCE  (FLAGtl,.F(l)) 

_ equivalence  ;w( i i >wave( i > i n . ( T( i i .tau'i  .1 n _ 

NSP= 144 

_ READ  (5,101)  t  W  ( K. )  ,  F  (  K  )  ,  ( COEF (  I  ,  K  )  >l  =  l>4)  iK-l.NSP) _ 

101  FORMAT  (F6.3tAl.3X.4E12. 3} 

_ DO  20  K= 1 >NSP _ 

FLAG! K)=. FALSE. 

_ If-LEXlUNf  «.Q  L-ELA&JXL=jlIRU£« _ 

20  CONTINUE 

1Q0Q  REAP  (5,102)  H20PCM.X02KM. PRESS. TAUSC _ 

1C2  FORMAT  (4F5.0) 

_ I£_ LE8JL&S*  dQ.Q.-l-  PREsS^I.. _ 

IF  (TAUSC.EQ.O. )  TAUSC=1 . 

_ C02ACM  =  X02XM«32. _ _ 

L( 1)=1./H20PCM 

_ L  t  2 )  =  1 . /SORT ( H20PCM*PR£SS  ) _ 

L ( 3 ) =1 ./C02ACM 

_ L<4  1=1. /SORT  (C02ACM«PR£SS) _  _ _ 

DO  10  K=1,NSP 

. . —1 F  <f  LAG  (  K )  ),.G0__TQ_14 _ _ 

DO  12  1=1.4 

12  A(  I )  =  (  L  (  I  )  »COEF  (I  ,K )  )  **2 _ _ 

XF  =  1 ,  /  SORT  (  A  ( 1 )  +  A  (  2  )  )  +1  •  /  SORT  (  A  (  3  )  -f  A (  4  )  ) 

_ GO  TO  11 _ _  _ 

14  XF=PRESS/COEF(3.K)/L(3) 

11  T(X!=TAUSC«EXP;-XF)«100.  _ _ _ _ _  _ 

IF  (XF.EQ.O.)  T(IC)=0. 

10  CONTINUE _ _ _ 

WRITE  (6,105) 

_ I.Q.L.f  QRMAl_Ll.HljJ_3A.iit3HJiQL--.£L[LAR._ABS.QfiEXLQN_Q£Ll[M£_INERAfi.ED  .lN  Alfi) 

WRITE  (6,106)  XO2KM.H20PCM 

106  FORMAT  ( 1H0, 13X. 13HPATH  LE NGTH  F5.2.4H  KM.I3X.F5. 2.11H  PR  CM  H 
120) 

_ WRITE  (6.108) _ 

108  FORMAT  llHO,  13X ,3 ( 16HWAVE-  TRANS-  ,6X)/ 

_ 11H  ,  1 3X .  3  ( 16HLENGTH  »*  MISSION, 6X)/ _ _ 

21H  » 13X » 3 ( 16HMICR0NS  PERCENT , 6X )/ 1H  ) 

_ WRITE  (-6.109)  (  (WAVE  (K,J) _ .TAUtK.Ji  .J*l  .31  .x=i  .<.») _ 

109  FORMAT  < ( 1H  ,13X,3(F6.3  , 2X , F6 . 1 »7X  )  ) ) 

_ WRITE  (6.119) _ _ _ 

no  FORMAT  (1H0.13X.63H**  WAVELENGTH  OF  LOCAL  MAXIMUM  OR  MINIMUM  OF  TR 

_ 1AN5MI SSION  CURVE) _ 

GO  TO  10C0  . 

_ JEND _ _ 

SDATA 

»**»»»  INSERT  ALL  144  ENTRIES  OF  TABLF  2  HERE  »»»»»» _ 

34.3  15.  1.  1. 
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»»»«»»« _ sample  outp u_t _ »♦♦»«* * 


MOLECULAR 

1  ABSORPTION 

_D£_THE  INFRARED. IN 

AIR 

PATH  LENGTH  15.00 

KM._  .... 

..  34,5.0__J>R  C.M..H20 

-  --  . 

WAVE- 

TRANS- 

WAVE- 

TRA  S- 

WAVE- 

TRANS- 

LENGTH** 

MISSION 

LENGTH** 

MISSION 

LENGTH**  MISSION 

MICRONS _ 

PERCENT 

...  MI  CRONS 

PERCENT 

MICRONS 

PERCENT 

3.000 

0.1 

3.535 

38.9 

..  4,3$Q 

0. 

3.010 

0.1 

3.540 

23.1 

4,450 

0. 

3.020 

0.1 

3.550 

25.0 

4.490 

0.4 

3.030 

0.2 

3.555 

23.3 

4.500 

0.0 

3.040 

0.0 

3.560 

43.4 

4.520 

0.0 

3.050 

0.0 

3.570 

28.6 

4.530 

0.0 

3.060 

0.0 

3.580 

39.8 

4.554 

3.5 

3.070 

0.0 

3.585 

22.5 

4.560 

2.7 

3.080 

0.1 

3.600 

47.7 

4.587 

12.7 

3.090 

0.0 

3.615 

31.2 

4.595 

9.6 

3.100 

0.0 

3.625 

61.4 

4.605 

15.5 

3.110 

0.0 

3.  635 

39.7 

4.615 

6.2 

3.120 

0.1 

3.645 

68*4 

4.630 

8.3 

3.130 

0.2 

3.650 

44.8 

4.  637 

6.0 

3.140 

0.1 

3.655 

49.8 

4.643 

6.8 

3.160 

2.4 

3.670 

9.3 

4.650 

5.9 

3.170 

1.8 

3.700 

69.4 

4.660 

7.9 

3.200 

0.0 

3.710 

37.0 

4.670 

1.5 

3.216 

0.0 

3.720 

63.0 

4.680 

12.9 

3.223 

0.0 

3.730 

31.1 

4.690 

8.7 

3.230 

0.1 

3.745 

57.5 

4.  700 

4.2 

3.240 

0.0 

3.755 

24.5 

4.715 

6.3 

3.250 

0.5 

3.  770 

47.  * 

4.  720 

2.7 

3.260 

0.1 

3.790 

36. 

4.730 

S.O 

3.270 

0.6 

3.  800 

70. 

4.  740 

6.2 

3.290 

0.0 

3.  8 10 

27.5 

4.745 

8.5 

3.300 

0.0 

3.  820 

67.1 

4.755 

3.3 

3.310 

0.0 

3.  830 

30.6 

4.760 

3.8 

3.327 

2.2 

3.840 

63.3 

4.780 

0.3 

3.333 

0.6 

3.  860 

37.9 

4.790 

2.3 

3.340 

0.7 

3.070 

53.8 

4.800 

0.1 

3.3  50 

0.3 

3.880 

40.0 

4.810 

0.0 

3.370 

4.0 

3.910 

71.2 

4.820 

1.0 

3.380 

4.3 

3.930 

45.1 

4.  830 

0.0 

3.400 

15.4 

3,940 

64.5 

4.840 

0.2 

3.402 

6.9 

3.950 

66.0 

4.855 

2.2 

3.420 

35.9 

3.960 

79.3 

4«860 

.1,5  .  ..  . 

3.427 

31.6 

4.000 

71.2 

4.870 

1.5 

3.433 

31.5 

4.050 

46.1 

4.890 

0.0 

3.440 

25.9 

4.  060 

43.1 

4.900 

0.0 

3.445 

37.0 

4.100 

26.0 

4.920 

1.2 

3.450 

26.3 

4.  140 

10.6 

4.930 

0.5 

3.46U 

41.7 

4.150 

7.6 

4.940 

0.0 

3.470 

37.5 

4.  160 

6.9 

4.955 

0. 

3.490 

40.9 

4.  180 

0.7 

4.965 

0.0 

3.500 

37.1 

4.  190 

0. 

4.,  975 

0.0 

3.514 

51.6 

4.350 

0. 

4.990 

0.0 

3.525 

41.0 

4.377 

0.0 

5.000 

0.0 

**  WAVELENGTH  OF  LOC AL  MAXIMUM  OR  MINIMUM*  OF  TRANSMISSION  CURVE 
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